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A mathematical model for the volcanic flow in planets is proposed. 

This mathematical model, which is one-dimensional, steady duct flow of a 
mixture of a gas and small solid particles (rock), has been analyzed. We 
apply this model to the lunar and the terrestrial volcanic flows under geo- 
metrically and dynamically similar conditions. Numerical results for the 
equilibrium two-phase flows of lunar and terrestrial volcanoes under simi- 
lar conditions are presented. The main results of our theoretical model 
are (1) the lunar crater is much larger than the corresponding terrestrial 
crater, (2) the exit velocity from the lunar volcanic flow may be higher 
than the lunar escape velocity but the exit velocity of terrestrial volcanic 
flow is much less than that of the lunar case. This result supports the 
hypothesis that Australian tektites came from the moon, as a stream of a 
mixture of rock and gas of extremely high speed and (3) the thermal effects 
on the lunar volcanic flow are much larger than those of the terrestrial case. 
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I. INTRODUCTION 



There is a real need for physical and mathematical models of geo- 
logic processes, especially in conjunction with planetological studies. 

The underlying reason is that we cannot observe planetological processes 
in action. This fact cuts us off from the traditional geologic method 
of explaining the structures of rocks in terms of processes which we see 
in operation. It is dangerous to try to explain planetary structures by 
comparing them with similar terrestrial structures, because we are extra- 
polating to very different physical conditions. For instance, lunar gravity 
is 1/6 of terrestrial; atmospheric pressure is 10 of the earth's; 
water content of the rock, undetectable; and oxygen fugacity (at^g^tic 
temperatures) 10 ^ of terrestrial. The extrapolation is so great that it 
is hard to be sure that a given lunar structure was produced by the same 
kind of process as an apparently similar terrestrial structure. In parti- 
cular, in view of the great differences in physical parameters, it is 
thoroughly unscientific to assume that the nearest analog of a terrestrial 
structure will be a lunar structure of about the same size. 

To some extent, it may be possible to replace the geological prin- 
ciple of uniformltarianism by a planetological principle of homology; we 
seek to explain planetary structures by comparison with terrestrial struc- 
tures, after allowance for scaling. There is a famous rule, called the 
pi-theorem (P.W. Bridgeman, 1922) which states that a physical law can be 
stated in terms of one or several non-dimensional parameters, each a pro- 
duct of a number of physical variables. If these numbers have the same 
value in two objects, one on planet A and the other on planet B, both pro- 
duced by the same process, then it is reasonable to try to explain what is 
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seen in B in terns of what is known about A. These two cases are in 
dynamical similarity. 

One of the most important problems of planetology is to understand 
the circular depressions which are the most conspicuous features of the 
surface of the moon, and the surfaces of most of the planets and satellites 
which have been examined. These depressions are all called craters; but 
it is important to realize that the use of one name does not mean that 
they are all the same kind of object on the moon or the planets, any more 
than on the earth, where some are due to impact and others to volcanism. 

J.F. McCauley (1967) has furnished a useful Indication of the truth 
of this observation, by drawing attention to the craters Maunder (Fig. 1) 
and Kopff (Fig. 2) in Mare Orientale. The two craters are of roughly the 
same size (Maunder 55 km, Kopff 40 km) but they are fundamentally different 
in structure. Kopff has a rim which is polygonal in plan and A-shaped in 
cross-section, while the wall of Maunder is circular in plan and much 
deeper on the inside than the outside; Kopff has a smooth surrounding sur- 
face, while that around Maunder is hummocky; Kopff is shallow while Maunder 
is deep. The point of McCauley's comparison is that both are on Mare 
Orientale, which has a peculiar, wavy . surface, clearly not degraded by 
any erosional process. The survival of this peculiar surface guarantees 
that neither crater can have been seriously modified by erosion; both must 
have been formed in about their observed shape. It seems to follow that 
the moon has at least two radically different way of forming craters. Since 
all impact craters tend to be similar, it follows that if Maunder is an 
impact crater, then Kopff is not. This example is important because it. 
indicates that the moon has large craters of volcanic origin. The explosion 
craters considered in this paper are probably different from Kopff in 

- ' : \ . 
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morphology. 

Similarly, W.S. Cameron and S.L. Padgett (1974) have drawn atten- 
tion to regional variations in the frequency of what they call ringdikes, 
i.e., craters more or less like Kopff. By a similar argument, C.R, Chapman 
et al. (1971) find that the most subtelescop ic craters on the maria are of 
internal origin, because they are ten times as numerous as the subtelescopie 
highland craters, whereas telescopic mare craters are one tenth as numerous 
as telescopic highland craters. 

In addition, the lunar origin of tektites, which seems especially 
probable for Australian and Javanese tektites, their flanges and ring- 
waves, seems to imply explosive lunar volcanism (O'Keefe, 1976); the argument 
rests in part on studies of australite ablation, especially those by D.R. 
Chapman and H.K. Larson (1963). The calculations have been thorougly 
checked by independent computations; they 6eem to rule out any non-lunar 
origin of tektites. The line of argument is supported by calculations 
of the size of the crater which would be needed for a terrestrial origin, 
and of the time needed to produce glass pieces of the necessary size, 
density and homogeneity, as well as of the probability of hitting composi- 
tions as unusual as those of some of the components of the tektite stream 
fields. 

On the other hand, in an elegant morphometric study, R.J. Pike (1 974) 
showed that most lunar craters resemble in shape terrestrial explosion pits, 
rather than terrestrial volcanoes. The exceptions are the maar-type craters 
and the tuff rings which are terrestrial, but are more like typical lunar 
craters than like the majority of terrestrial volcanoes. Pike rejects the 
the maars and tuff rings, because the maximum size of these objects is about 
5 km, while lunar craters range up to 250 to 300 km. Obviously the question 


of scale is crucial here 


It thus appears useful to calculate the physics of eruptive volcanism, 
with the earth and the moon as examples. (Martian eruptive volcanism will 
then represent an intermediate case.) 

In comparing lunar and terrestrial volcanism, we note that water, 
which is the principal driving gas in a terrestrial volcano, will be less 
important than hydrogen in a lunar volcano. This is a direct consequence of 
the low oxygen fugacity in a lunar lava. The equilibrium is governed by 
the usual relation 


(h 2 o) 

(H 2 ) (Oj )* 9 


( 1 ) 


where (H^Q) is the partial pressure of steam, etc. , and K is the equi- 
librium constant; from the JANAF tables (Stull and Prophet, 1971) we find, 

-2 

at 1500*C log^Q K - 3.225 if pressure is in N m (Newtons per square meter). 

The fugacity of oxygen in lunar lava at this temperature is known to be about 
-13 -2 

10 N m ; thus is the is the partial pressure of oxygen in contact with 
the lava. Substituting this in Eq. (1), we find that the ratio of the partial 
pressures of water and hydrogen is 0.17. Hence in our mathematical model, 
we consider steam as the gas in terrestru l volcano while hydrogen as the gas 
in lunar volcano. 


II. MATHEMATICAL MODEL OF VOLCANIC FLOW 

To explain the mechanism of volcanic eruption, whether lunar or ter- 
restrial, we may consider the volcano as a two-phase system consisting of the 
solid particles and a gas. We analyze the one-dimensional flow due to the 
eruption of such a two-phase system. From the field observation (G. Fielder, 
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1971) and some preliminary theoretical calculations [T.R. McGetchin and 
G.W. Ullrich, 1973], the general volcanic flow field is sketched in Fig. 3. 

The general flow field may be divided into two regions: the lower 

part, region I, starts from a great depth in the planet; region I has a 
height of the order of 100 km in terrestrial case. At the base, due 
to chemical reaction, the rock is melting at a temperature of bhe order of 
100 *C. The original pressure is very large [T.R. McGetchen and G.W. Ullrich, 
1973]. As a result, a mixture of molten rock and gas flows upward in a 
duct of almost constant cross-sectional area. In this region, since the 
pressure is very high, the mixture of the molten rock and the gas behaves 
like an incompressible fluid. From fluid dynamic point of view, the flow 
problem is very simple. The flow field should not deviate much from the 
simple analysis by T.R. McGetchftn and G.W. Ullrich [1973]. We shall not 
discuss it any more. 

When the flow is near the surface of the planet, with x - H Q , 
where x is the vertical coordinate, the area of the duct Increases and the 
mixture expands as It flows upward. We call this region as region II, in 
which the flow field is very complicated. In this paper, we study the es- 
sentlal features in region II of lunar and terrestrial volcanic flows. We 
are especially interested in the similar flow patterns of a terrestrial 
volcanic flow with a corresponding lunar volcanic flow. 

III. FUNDAMENTAL EQUATIONS OF VOLCANIC FLOWS AND THE IMPORTANT NON-DIMENSIONAL 
PARAMETERS. 


Since the flow of the mixture comes a long way from the great depth 
of the planet to the initial section x ■ - H q of region II, it should reach 
the equilibrium condition such that the initial velocity of the rock and that 


of the gas are approximately equal and that the temperatures of the rock 
and the gas are also equal. It is known as equilibrium flow. (S.I. Pai, 
1977). In order to get some essential features of the volcanic flow in 
region II , we assume that the flow in region II remains as equilibrium flow 
first. We shall investigate the non-equilibrium effects on the flow field 
in region II in a later paper. 

For equilibrium 'flow of the two-phase system in region II, we 
consider the following variables: 

velocity of the mixture u 
temperature of the mixture T 
pressure of the mixture p , and 
density of the mixture p 

The density of the rock P r is assumed to be constant and the 

density of the mixture depends on the density of the gas p and the volume 

S 

fraction of the rock in the mixture Z . 

The values of the variables at the initial section x * - H of 

o 

cross-sectional area A q are given from the calculation of region I as 
follows: 


p - P o , u - u Q , T - T q , P - P 0 • (2) 

Since we are in- -rested in similar flows in the lunar and the terres- 
trial case, the relations between the initial values of these 
two cases may be obtained by similarity considerations as we shall show 
in section IV. 

We define the non-dimensional quantities (with a bar) as follows; 



where L is the characteristic length of the flow field of region II 
which may be taken as H q • 

The sound speed of the gas a is given by the following formula 


a « 



(4) 


where y is the ratio of specific heats c p/° v » \ t * le universal 

gas constant and m is the molecular weight of the gas. 

The non-dimensional equations for one-dimensional steady equi- 
librium flows of a two phase mixture of solid particles and a gas (s.I. Pai 
1977) are as follows: 

(i) Equation of state 



where subscript o refers to values at initial section x 
(ii) Equation of continuity 


(5) 


H 


o 


p u A ■ 1 

(iii) Equation of motion 



( 6 ) 


( 7 ) 


where the sound speed of the mixture a^ is given by the formula; 


where 


2 


™ - V 

(1 “ Z) 2 




( 8 ) 


( 9 ) 


is the mass concentration of the solid particles in the mixture and P 
is the effective ratio of specific heats of the mixture, i.e., 



1 + r\& 


( 10 ) 


where n-l-h/k , <5 ■ c /c ,c is the specific heat of the solid 

p p s v s 

particles, c v is the specific heat of the gas at constant volume and c^ 

is the specific heat of the gas at constant pressure. 

The initial Mach number M is defined as follows; 

o 


M 

o 


u 

o 

a 

o 


and the Initial Froude number Fr is defined as follows; 

o 

2 


u 



(ID 


( 12 ) 


where g is the gravitational acceleration of the planet. The friction 
coefficient of the duct is c f and 5 - D/L where D is the hydraulic 
diameter, A(x) is the cross-sectional area of the duct in region II. 


(iv) The energy equation: 




lb <1-Z) 2 + b 4 -^1 + b <X-Z> 2 J 

uA o 

- I f- (i a+1 - - 

n-o D 

av 

U. V^-lJa^V C 0 p g0 r 0 (Y.-1) r o J 2 C Fr 0 


a - k ) (y - «r 

p 


(y - 1) r 


b . il i- A 
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is the initial density ratio and D &v is the average hydraulic diameter 
between two steps and u^ v is the average velocity between two steps 
(Y. Hsu, 1976). 

Prom the non-dimensional equations (5) , (6) , (7) and (13) , we 

have four Important non-dimensional parameters: (9); M q (11) ; 

Fr (12) and G (15), which characterize the volcanic flow. The si£- 
o 

nificance of these parameters is as follows. 

The mass concentration k indicates the mass fraction of the 

P 

solid particles in the mixture, which is important in the determination 
of the velocity of the mixture. 
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The initial density ratio G indicates the relative importance 
of the gas and the solid particles in the mixture. 

The Mach number M q is a measure of the compressibility effect 
of the mixture due to high speed. 

The initial Froude number Fr Q indicates the effect of gravi- 
tational force on the flow field. 

IV. SIMILARITY LAWS OF LUNAR AND TERRESTRIAL VOLCANIC FLOWS 

We would like to compare a terrestrial volcanic flow with a 
similar lunar volcanic flow. For the similitude of these two volcanic 
flows, we should have both geometrical similitude and dynamic similitude 
(L.I. Sedov, 1959). 

By geometrical similitude, we mean that there is a linear scaling 

of the region II of the terrestrial and the lunar volcanoes. We should 

compare these volcanic flow fields with the same non-dimensional initial 

area and same slope of the non-dimensional variation of the cross- 

sectional area in region II, i.e., the same dA/dx . These two volcanoes 

are geometrically similar. It should be noted that the ratio of the exit 

area A to the initial area A for these two volcanoes are not the 
s o 

same because we have to satisfy the conditions of dynamic similarity. 

By dynamic similitude, we mean that the most important non- 

dimensional parameters of the fundamental equations of the volcanic flows 

for these two cases whould be equal. In other words, for dynamic similitude, 

the non-dimensional parameters: M , G , Fr and k for the similar 

o o p 

terrestrial and lunar volcanic flows should be equal. From these condi- 
tions, we may determine the corresponding size and initial values of these 
two similar volcanic flows as follows. 




We assume that in ths tarrastrlal volcanic flow, tha gas is 
steam vhlle in the lunar volcanic flow the gas is hydrogen as we have 
discussed in section I. We have the following values of various charac- 
teristic quantities in two geometrically and dynamically similar terres- 
trial and lunar volcanic flows: 

(i) For geometrical similitude, we have in region II 


A - A " 1 , 
om oe 


(& -(f) 


where subscript o refers to the value at the initial section x ■ - H 
or x ■ - 1 , ^ subscript m refers to the v*^.ue for lunar volcanic flow 
while subscript e refers to the value for terrestrial volcanic flow. V: 
compare two similar volcanic flows that satisfy Eq. (16). 

(ii) We assume that the initial temperature T 0 is the same 
in these two similar tarrastrlal and lunar volcanic flows, i.e.. 


X - T 1,000 *C 
om oe 


. (iii) In the terrestrial volcanic flow, the principal gas is steam 
while that in the lunar volcano is hydrogen. Hence at T^" 1000*C, we have 

c 

m - 18.016 ; c - 0.593 cal/gm *K ; Y ■ “ 1.23 (for steam) (18a) 

a P« • c ve 

c 

n - 2.016 : c - 3.721 cal/gm *K ; y - - 1.36 (for hydrogen). (18b) 

in pm m c 


The ratio of the initial sound speed in these two cases is 


- 3.126 




M i 1 1! i iMnijKiip'Miii 1 !!' |] 




(lv) By dynamical similarity, the Initial Mach number of these two 
similar flows should be the same, i.e.. 


om a 


Hence 


3.126 u 


For dynamical similarity, the initial value of velocity in the lunar case 
should be 3.126 times that of the terrestrial case. 

(v) By dynamical similarity, the density ratio 0^ of these two 
similar flows should be the same, i.e., 


om p 


oe p. 


Hence 


P “ P 
gom goe 


because we assume that the rock densities p_ in these two cases are the 

■r 

same. Since the initial density of the gas is related to the initial 
pressure of the mixture by the following formula 


go Ru 


m o 


we have the relation of the two initial pressures for those two similar 
volcanic flows as follows. 


e _ 

p - — p & 9 p 

om m oe oe 

m 


( 25 ) 
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If we assume that the pressures in the reservoirs are the same and because 

the smaller gravitational forces, the drop in pressure as the flow ascends 

must be smaller in the lunar case. Hence it is possible that p is larger 

om 

than p 

r oe 

(vi) By dynamical similarity, the initial Froude numbers Fr Q of 
these two similar flows should be the same, i.e., 

2 2 

u u 

TV . ...pm - m ff r . --PP— 

om g L * r oe g L . (26) 

mm e e 


Hence for dynamic similarity, the ratio of the two characteristic lengths 
of these two similar volcanic flows are 


L 

u 

m m 

om 

L 

2 

e 

u 


oe 



(27) 


We shall take L ■ H 

o 

(vii) By dynamic similarity, the should be the same in these 

two similar volcanic flows. In our number calculations, we calculate the 

flow field corresponding to various values of k^ and then compare the 

terrestrial and lunar volcanic flow with the same k . 

P 

For simplicity, in our numerical calculations, we shall assume that 


— ■ K « constant. (28) 

dx 

We calculate the cases for K ■ 1 , 3, and 5 and the corresponding ratio 
of the exit area A g to the ir^ial area A q are 

K - 1 , 3 , 5 

<We -2,4, 6 

(V A 0 > m " 61 , 181 , 301 . 
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Since ve are mainly interested in the flow variables at the exit of the 

volcano, the ratio A /A is more important than the slope dA/dx which 

s o 

determines the flow inside the duct of the volcano. It is interesting to 

notice that the exit area, the crater area of the lunar case, is about 

30 to 50 times that of the terrestrial crater if the initial areas of the 

two volcanoes A are the same, 
o 

V. NUiwRICAL SOLUTIONS FOR SIMILAR LUNAR AND TERRESTRIAL EQUILIBRIUM 
VOLCANIC FLOWS 


We shall show later that the viscous effects in region II are 
negligible (see Fig. 14). Hence in our numerical calculations, we con- 
sider the equilibrium inviscid flow only so that the friction terms in 
Eqs. (7) and (13) are neglected. 

The initial conditions in our numerical calculations are: 

T Q - 1,000°C , M q - 1.0 , c s - 0.3 cal/gm°K , p r - 3.30 gr/cc 

p ■ 100 atmospheres, L^ “ H * 1 km 
r oe r ’ e oe 

p ^ 900 atmospheres, L ■ H - 60 km 
r om r m om 

and the values for the gases are given in Eq. (18). 

We have calculated two cases. 

(1) Isothermal case 

In this case, we assume that the temperature T of the mixture 
remains at its initial temperature T q . This is physically possible 
because when the mixture moves upward, some of the molten rock may be 
solidified and release sufficient heat so that the resultant temperature 
of the mixture remains unchanged during the region II. Furthermore, from 


our previous experience in ash flow calculations, we find that the isother- 
mal approximation is a good approximation for the terrestrial case (?ai, et al. 
1972). We would like to analyze this simpler case first. 

(2) Adiabatic case 

In this case, we assume that there is no heat transfer between the 
* 

flow and the wall of the duct, nor viscous dissipation. The temperature 
of the mixture is not a Constant in region II but decreases according to 
the isentropic process as it ascends. 

Detailed numerical results are given in the reference by Y. Hsu 
(1976). Here we just give some typical results in Figs. 4 to 14 to show 
the similarities and the essential differences between similar lunar and 
terrestrial volcanic flows. 

In Fig. 4, we compare the actual velocity in the isothermal similar 

lunar and terrestrial volcanic flows with the mass concentration k as 

P 

a parameter and the corresponding area ratios (A fl /A o ) e ■ 4 and (A g /A o ) o *i81, 

The lunar crater area is about 45 times that of the similar terrestrial' 
one in this case. The velocity in all cases increases as k^ decreases. 

The velocity in the lunar case is much larger than the corresponding value 

in the terrestrial case. When k is equal to or less than 0.8, the exit 

P 

velocity of the lunar volcanic flow is larger than the escape velocity on 
the moon which is 2,372 m/sec. 

Fig. 5 shows the effect of the area ratio (A /A ) on the velocity 

so 

distribution in the Isothermal case. When A /A increases, the velocity 

s o 

increases too, but the effect of (A /A ) on the velocity is smaller than 

so 

that due to k . In the lunar case, the velocity increases rapidly in the 
P 

first part duct and remains almost constant for the portion of the 

duct. 

Fig. 6 shows the effect of k^ and (A^/A^) on the pressure 


distributions for the isothermal terrestrial volcanic flow. The influ- 


ence of k on the pressure distribution is small, particularly when 
P . 

(A /A ) is large, 
s o ° 

Fig 7 shows the effect of k and (A /A ) on the pressure 

p 6 O 

distributions for the Isothermal lunar volcanic flow. The Influence of 

kp on the pressure distributions is negligible. 

Figs. 8 to 13 show the results of adiabatic volcanic flows. In Fig. 

8, the velocity distributions of the similar lunar and terrestrial volcanic 

flows with corresponding area ratio (A /A ) "4, and (A /A ) * 181 

soe s o m 

are given with as a parameter. In general, these .velocity distri- 
butions are similar to those of isothermal case (Fig. 4) but the velocity 
is less than the corresponding value of the Isothermal case. The initial 

rate of increase of velocity of the lunar case is larger than that of the 
isothermal case. For small value of k . , the velocity of the lunar vol- 

canic flow will reach its maximum possible value at a short distance 
above the initial section. From that point up, the lunar volcanic flow 
will be a jet of a given width less than the width of the duct moving with 
the constant velocity of the value of maximum possible value. 

Fig. 9 shows the temperature distributions in the terrestrial 
adiabatic volcanic flow. The temperature drops very fast initially, especi- 
ally when kp is small . In the major portion of the duct, however, the tem- 
perature is almost constant but at a smaller value than- the initial temperature. 

Fig. 10 shows the effects of k and (A /A 5 on the tempera- 

p so 

ture distributions of the adiabatic terrestrial volcanic flow. For large 
kp ( > 0.8) , the Isothermal approximation is a good one for terrestrial 
volcanic flow. 

Fig. 11 shows the effects of k and (A /A ) on the temperature 

" p s o 


i.3 

distribution® of th« lunar adiabatic volcanic flow. The thermal affect 
on the lunar volcanic flow is larger than that of terrestrial flow. The 
Isothermal approximation is not good for lunar volcanic flow. The tempe- 
rature of the lunar adiabatic flow drops very fast initially, especially 

when k is small. When k is less than 0.5, the condition of m a x im u m 
P P 

possible velocity and zero temperature will be reached at a short distance 
from the initial section. 

Fig. 12 shows the effects of k and (A /A ) on the pressure 

p s o 

distributions for the adiabatic terrestrial volcanic flows while Fig. 13 

shows the corresponding pressure distributions of similar lunar adiabatic 

0 

volcanic flow. The pressure drops initially very fast when k^ is small 
but remains almost constant in the major portion of the duct. The pressure 
in the lunar case is much smaller than the corresponding value in the lunar 
case even though the initial pressure of the lunar case is larger. 

Finally, Fig 14 shows the effect of friction on the velocity 
distribution on the isothermal terrestrial volcanic flow. The effect of 
friction is negligible. 


VI. SUMMARY AND CONCLUSIONS 

From our theoretical study and numerical results, the following 
conclusions may be drawn. 

(1) The flow fields of similar terrestrial and lunar volcanoes 
have been calculated by the two-phase flow theory for a mixture of a gas 
and solid particles, considering both isothermal and adiabatic equilibrium 
cases. 

(2) The important non-dimensional parameters which characterize 
the volcanic flow are (1) the Mach number (ii) the Froude number, (Hi) the 



of the rock in the volcanic flow. 

(3) By dynamical similarity, the above four non-dimensional para- 
meters should be the seme. If hydrogen powers lunar volcanoes, this suggests 
that terrestrial volcanoes with steam as the gas are to be compared (other 

things being equal) with lunar volcanoes 60 times larger in area. Hence the morpho- 
logical similarity found by Fike (1974) between terrestrial maars and tuff 
rings, on the one hand, and the commonest kind of large lunar crater on the 
other is not necessarily invalid, despite the fact that the lunar craters 
are SO times larger than the maars and tuff rings. 

(4) In the terrestrial cases, the exit velocity from the volcano 
is supersonic, but its value is much smaller than the terrestrial escape 
velocity. 

(5) In the lunar cases, the exit velocity from the volcano may 

be higher than the lunar escape velocity when the mass concentration of the 

rock is less than k n - 0.8. 

P 

(6) In the adiabatic case, the lunar volcanic flow may reach its 
maximum possible velocity before the exit and the volcanic flow will then 
leave the volcano as a jet stream with velocity greater than the lunar es- 
cape velocity. Thus the flow field agrees with the assumed flow pattern 

of a stream of tektites coming from the moon (D.R. Chapman and H.K. Larson, 

1963). 

(7) When k is large and (A /A ) is small, the temperature in 

p so 

the terrestrial cases is quite close to its initial value. When k is 

P 

small, the initial drop of the temperature in both the terrestrial and the 
lunar cases is large; in the major portion of the duct, the temperature 
is almost constant but is less than the initial temperature* 



(8) The variation of both the velocity and the pressure in 
the terrestrial cases is gradual. For the lunar cases, the pressure drops 
very rapidly to a very low value initially and then remains almost con- 
stant for the major portion of the duct. The variation of the pressure 
with k n is always negligibly small. 
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VIII. CAPTIONS AND FIGURES 

Fig. 1. The lunarcrater Maunder, 55 km in diameter, on Mare Orientale. Note 

the hummocky region behind it, and the very gentle slope of the outside 
wall* The slope of the inner wall has been diminished by slumping. A 
smaller crater, Hohmann, near the bottom of the figure, is analogous 
to Kopff . 

Fig. 2. The lunar crater Kopff, 40 km in diameter, in Mare Orientale. Note 
the steep outer wall, and the smooth region around it. Kopff is ob- 
viously a different kind of crater from Maunder, yet it is superposed 
on the characteristic Mare Orientale Terrain, which has not been re- 
worked by cratering. It is thus not likely that Kopff is the worn- 
down remnant of a crater like Maunder. 

Fig. 3 Mathematical model for similar lunar and terrestrial volcanic flows. 
Comparison of the actual velocity of i sothermal similar lunar and 
terrestrial volcanic flows along the duct with the k^ as a parameter, 


Fig. 4. 
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for the corresponding area ratio (A /A ) • 4 and (A /A ) 

a o a a ora 

- 181 . 

Fig. 5. Comparison of tha actual velocity distributions of isothermal 

similar lunar and terrestrial volcanic flows along the duct with 

the corresponding area ratio as a parameter for a fixed k * 0.5. 

• P 

Fig. 6* Sf facts of k and area ratio po the non-dimensional pressure 

P 

distributions of Isothermal terrestrial volcanic flow along the 
duct. 

Fig. 7. Effects of k^ and area ratio on the non-dimensional pressure 
distributions of Isothermal lur : volcanic flow along the duct. 

Fig. 8. Comparison of the actual velocity of adiabatic similar lunar and 

terrestrial volcanic flows along the duct with the k p as a para- 
meter for the corresponding area ratio (A g /A o ) e ■ 4 and (A s /A Q ) a 
- 181. 

Fig. 9. Temperature distribution of the adiabatic terrestrial volcanic flow 
along the duct with as a parameter and (A g /A Q / fi " 4. 

Fig. 10. Effect of k p and area ratio on the temperature distribtuion of 
terrestrial volcanic flow along the duct. 

Fig. 11. Effects of kp and area ratio on the temperature distribution of 
lunar volcanic flow along the duct. 

Fig. 12 Effect of k and area ratio on the non-dimensional pressure 

P 

distribution of adiabatic terrestrial volcanic flow along the duct. 
Fig. 13. Effects of k p and area ratio on the non-dimensional pressure dis- 
tribution of adiabatic lunar volcanic flow along the duct. 

Fig. 14. Effect of friction on isothermal terrestrial volcanic flow. 




» 1 



R FPRODUCIBILITy OF THE 
PAGE 18 POOR 










I0 4 

9.47 


tyo. 


Kp=0.5 


K p *0.8 
MOON Ves 


Kd s 0.95 



Kp s 0.5 


K p *0.8 


K„=0.95 


_ MOON, A S /A 0 =I8I 

EARTH, A s /A 0 =4 

1 1 i ' i 

0 -0.2 -0.4 -0.6 -0.8 


Fia. 4 











EARTH 


VQ-'V// jj 

A s /A 0 =2 K d =0.5^^ // 

Kp-O.^ j 

: p=°. 95 _^^ / / 

S^^/A 0 =4 

Ko=0.95 *-e?**^ 


'^0=0.95-0.1, A,/Ao=6 






ISOTHERMAL 



K p =0.95 

ISOTHERMAL 

















VISCOUS 

EARTH, A s /A 0 =4 







